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Abstract 


This analysis paper presents previously unknown properties of some special cases of the 
Wright function whose consideration is necessitated by our work on probability theory and 
the theory of stochastic processes. Specifically, we establish new asymptotic properties of 
the particular Wright function 


OO 


r(fc + pn) x’^ 
r(pn) n! 



(1*1 < oo) 


when the parameter p G (—1, 0)U(0, oo) and the argument x is real. In the probability theory 
applications, which are focused on studies of the Poisson-Tweedie mixtures, the parameter 
fc is a non-negative integer. Several representations involving well-known special functions 
are given for certain particular values of p. The asymptotics of i’Li(p, k; p, 0; x) are obtained 
under numerous assumptions on the behavior of the arguments k and x when the parameter 
p is both positive and negative. 

We also provide some integral representations and structural properties involving the 
‘reduced’ Wright function o^'i(—; P, Oi *) with p € (—1,0) U (0, oo), which might be useful 
for the derivation of new properties of members of the power-variance family of distributions. 
Some of these imply a reflection principle that connects the functions o^i(—;ip, 0; •) and 
certain Bessel functions. 

Several asymptotic relationships for both particular cases of this function are also given. 
A few of these follow under additional constraints from probability theory results which, 
although previously available, were unknown to analysts. 
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pansions, Stokes lines, reflection principle, multiplication property, Airy function, Bessel 
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1. Introduction 


In this analysis paper we deal with various asymptotic properties of a particular special function 
known as the Wright function that arises in analytic probability theory. This article has its 
outgrowth in the specialist works of [5] and [29] and is closely related to [14, Sections 3, 4]. 
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We consider the subclass of Wright functions defined by 


z) = '^ 


n—0 


r(an + a) z" 
T{j3n + b) nl 


( 1 . 1 ) 


where the parameters a and /3 are real and positive and a and b are in general arbitrary complex 
numbers. We also assume throughout that the a and a are subject to the restriction 


an + a ^ 0, —1, —2 ,... (n = 0,1, 2,...) 


so that no gamma function in the numerator in ( 1 . 1 ) is singular; however, if the gamma function 
in the denominator is singular at the same points as the numerator, then a limiting procedure 
may be applied. In the special case a = /3 = 1, the function i'I'i(z) reduces to the confluent 
hypergeometric function (r(a)/r(&)) iFi{a]b; z); see, for example, [20, p. 40]. When a = I, the 
function i4'i(z) becomes the generalized Mittag-Leffler function considered in [19]. 

We introduce the parameters associated with (1.1) given by 

K = 1 + (3 — a, h = , d = a — b. (1-2) 

If it is supposed that a and (3 are such that k > 0 then i'I'i(z) is uniformly and absolutely 
convergent for all finite z G C. If r = 0, the sum in (l.I) has a finite radius of convergence equal 
to h~^, whereas for k < 0 the sum is divergent for all nonzero values of z. The parameter n is 
found to play a critical role in the asymptotic theory of 1411 ( 2 :) by determining the sectors in the 
complex z-plane in which its behavior is either exponentially large or algebraic in character as 
|z| —>■ 00 . 

The consideration of the entire class (1.1) of Wright functions is interesting in its own right 
and has already been undertaken in [14, 26]. In contrast, in this paper we concentrate on the 
derivation of subtle properties of a few special cases of this function. This is partly motivated by 
the fact that the function in (1.1) with a = j3 = p,b = 0 and a = k, where k denotes a non-negative 
integer, has numerous applications in probability theory, the theory of generalized linear models 
and stochastic models of actuarial science. This is because it turns out that the probability 
function of a generic member of the specific three-parameter family of univariate probability 
distributions on non-negative integers, which are hereinafter referred to as the Poisson-Tweedie 
mixtures, is expressed in terms of the function i4'i(z). We refer to [4, 6 , 7, 9, 13, 16, 17, 23] for 
more detail on this important family of discrete probability distributions. The simplest cases, 
which correspond to the values P = — 5 and p = 1, are considered in [23] and [16], respectively. 
The observation that the probability function of a generic member of this family is expressed in 
terms of i'I'i( 2 :) appears to be new and will be pursued in depth in our forthcoming paper [17]. In 
contrast, [13, Section 7.5] and [ 6 , §11.1.2] employ a double recursion for the specific probabilities 
associated with Poisson-Tweedie mixtures. 

In this paper, we consider the particular case of the function in (1.1) given by 


i'^iip,k;p,S;x) = ^ 

n—0 


r(A: -I- pn) x" 
r(pn-|-(5) nl 


(fc = 0,l,2,...), 


(1.3) 


where S G C, p £ (—1,0) U (0,oo) and throughout the paper x is a real variable. Although k is 
a non-negative integer in the probability application, we shall occasionally relax this constraint 
and consider situations where k > 0 is not necessarily an integer; see, for example, (1-5). 

In a similar manner to [18, Section 2.3], we also define the complex-valued ‘reduced’ Wright 
function (j)(p, S, z) (also known as a generalized Bessel function) of argument z £ C, which is 
indexed by the real-valued parameter p £ (—1,0) U (0,oo) and the complex-valued parameter 
5 £ C hy the following convergent series: 


]4'i(—;p,(I;z) = ^(p,5;z) := ^ 


n—0 


nl T{pn + 6) 


(1.4) 


In this paper, we shall set the parameter d = 0 in (1.3) and (1.4). In this case, the first term of 
the series corresponding to n = 0 on the right-hand sides of (1.3) (when k 0) and (1.4) is zero; 
compare to formula (5.1). 
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The derivation of the new results of the present paper, which concern numerous properties of 
members of the classes of the special functions (1.3) and (1.4) (with 6 = 0), is necessitated by 
our studies of various approximations for some Poisson-Tweedie mixtures and members of the 
power-variance family of probability distributions. In view of [24, 25], numerous representatives 
of this power-variance family are closely related to the ‘reduced’ Wright function (1.4) with 5 = 0. 

At this point, we observe the following important relation for a; > 0 between the corresponding 
representatives of the two classes (1.3) and (1.4) of special functions introduced above. This is 
given by 

pOO 

p,0;ex) = / <p{p,0;exT^)dT, e = signp (1-5) 

Jo 

valid for k > 0 (not necessarily an integer) when p > 0 and for non-negative integer values of k 
when p G (—1, 0). The proof of (1.5) is given in Appendix A. 

The structural relationship (1.5), which connects the Wright function with the corre¬ 
sponding ‘reduced’ Wright function (f>, provides the analytical explanation of the fact that the 
so-called Poisson-Tweedie mixtures are derived by employing the probability law of a member of 
the power-variance family as the mixing measure. See also [23, Eq. (16)] and [16, Eq. (3.9)] which 
pertain to the values p = — 5 and 1 , respectively. 

The structure of the paper is as follows. In Section 2, some special function and polynomial 
representations of i4ii(p, k; p, 0; x) are given. In Sections 3 and 4 the asymptotics of this function 
are considered for (i) x —>■ ± 00 , k finite (ii) k —>■ -l-oo, x finite when p > 0 and p G (— 1 , 0 ), (iii) 
k —>■ -1-00, X —>■ -1-00 when p > 0, (iv) k —>■ -l-oo, x —>■ —00 when p G (—1,0) and (v) p —>■ -boo when 
fc is a positive integer and x is finite. The asymptotics of the ‘reduced’ Wright function ^(p, 0;a;) 
for X —>■ ±00 are summarized in Section 5. The concluding Section 6 lists various representations 
and properties of (j){p,0]x) when p G (—1,0). The Appendices A, B and C contain a proof of 
(1.5), two integral representations for i'i>i{p,k-, p,0;x) stated in Section 2 as formulas (2.5)-(2.6) 
of Theorem 1 and two specihc representations of the ‘reduced’ Wright function in terms of Bessel 
and Whittaker functions. 

To conclude the Introduction, we point out that the majority of the results of this paper 
provide an analytic foundation to, or constitute extensions of the analysis counterparts of, the 
corresponding assertions of probability theory and the theory of stochastic processes. A few of 
the analytical results presented only have probabilistic proofs; in these cases, as for many of the 
other results presented herein, a numerical verification has been carried out. 


2. Special representations of iTi(p, A:; p, 0; x) 


For certain values of the parameter p it is possible to give representations of i4ii(p, fc; p, 0; x) in 
terms of well-known special functions. The simplest of these pertain to the values p = 1 and p = 
— |, —i, —|. See Sections 2.1-2.2 and also Sections 6 .1-6.3 for the case of the ‘reduced’ Wright 
function (j){p,0;x) when p = 1 and p = —j, —|, i, i, |. We remark that the representations in 
Sections 2.1 and 2.3 are given for real x, but are valid more generally for complex values of x. 
2.1 Special function representations 
When p = 1 and fc > 0, we have 


i4'i(l,fc; l,0;x) 


^ r(fc -(- n) x” 

^ r(n)n! 

n—1 ^ ' 

xk\ iFi{k + 1 ; 2 ; x), 


00 


E 

n—O 


r(fc -b n -b 1 ) x" 

r(n -b 2 ) n! 


( 2 . 1 ) 


where iFi denotes the confluent hypergeometric fnnction. In view of (2.1), numerous analytical 
results of our paper [16] constitute special cases of related assertions of the present article. 

Another case where a special function representation is possible corresponds to p = In 
this case, we have 


i^'i(-i,fc;-i, 0 ;±x) 


^ T{k - in) (±x)” 

•^T{k — n) x^" ^r(A: —n—^) x^”+^ 

r(-n) ( 2 n)! ^ r(-n-i) ( 2 n-bl)! 


( 2 . 2 ) 


El ±S2. 
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Then, for positive integer fc, some straightforward manipulation yields 


Si 


n—k 


r(l-fc + n)r(n+i) 


m=0 


(la;)2m+2fe 


TO!r(fc + i + m) 


and 


1 °° 

n 

n=0 




n\ 




T(fc-n-i)(ixf"+i = (-)'= 

n =0 


{\x) 


2n+l 


n!r(| — A: + n) 


where Iv{x) denotes the modified Bessel function of the first kind. Consequently, for non-negative 
integer /c, we have the representations 


i4'i(-i,fc;-i,0;a:) = (2.3) 

and 

i4/i(-i, fc; -i 0; -x) = (2.4) 

where we have employed the well-known relation (x) = (x) -I- ( 2 /tt) sin( 7 ri^) K^, (x) connecting 

I±v{x) with the modified Bessel function of the second kind K^{x) (also known as the Macdonald 
function); see, for example, [11, Eq. (10.27.2)]. 

It is worth mentioning that when k > 0, but excluding half-integer values (for which the 
left-hand side of ( 2 . 2 ) is not defined), then Si = 0 and we easily find that for a; > 0 


i'I'i(— 5 , fc; — i, 0 : ±a;) = ± 


cos nk ^ ^ 


.(x) 


(fc>0; i,l,|,2,...). 


2.2 Integral representations when p = —1/3 and p = —2/3 

The Wright function i'I'i(p, fc; p, 0; —x) when p = —^ and p = —| is also related to special 
functions, but in these cases it is expressed as an infinite integral involving the Ki, Bessel and W^.,^ 
Whittaker functions, respectively. Originally, we obtained these results by using probabilistic 
arguments. In this paper, we provide an analytical proof. These are given in the following 
theorem. 


Theorem 1. For x > 0 and positive integer k, we have 


and 


i4'i(-i,A:;-i,0;-x) 


Stt 


pOD 

Jo ^ 


2x3/2 

3v^ 


dt 


(2.5) 


i4'i(-|,fc;-|,0;-x) 


e 3 exp 



Wi 1 
2 ’ 6 



( 2 . 6 ) 


The proofs of (2.5) and (2.6) are given in Appendix B. We should point out that for p = 
1,—j,— 5 ,—|, the results in (2.1), (2.4), (2.5) and (2.6) can all be derived from (1.5); see also 
Section 6.1. 


2.3 Polynomial representation for p ^ 0 

It is important that, for positive integer k and p ^ 0, the Wright function i'I'i(p,/c; p, 0; x) can 
be expressed as (px)^e® multiplied by a polynomial in x~^ of degree k—1. This is established in 
( 2 . 8 ), which is the main result of this subsection. 

To see this, we first note that the quotient of gamma functions 

k 

pn(l -f pn)... (fc - 1 -f pn) = {-f ^ si'"^(-pn)’” 

m—0 

k k 

E n(n - 1)... (n - r + I), 

m=0 r=0 


r{k + pn) 
T{pn) 
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where and are the Stirling numbers of the first and second kinds, respectively (see, for 
example, [11, p. 626]); some Stirling numbers are given at the beginning of Section 3. Then 


i5'i(p,A:;p,0;x) = ^ 


r(fc + pn) X 
—^ r(pn) n! 




m—O r—0 


(n — r)! 


k m 




m—O r=0 


where the sum over n has been evaluated as . 

If we put k — r = n and k — m = i, where 0 < £ < k and i < n < k, then the above double 
sum can be rearranged as 

k k k—1 n 

1=0 n=i n=0 e=o 


where, since Sm^ = Sm^ = (5o,m, the inner sum vanishes when n = k. If we define the polynomial 
hm(u) of degree m by 

m n 

h,„(u) := ^ := (2.7) 

n=0 £=0 

then, for all x (both positive and negative) and p ^ 0, we can write 

i4'i(p, fc;p,0;x) = (px)''e''hfc_i((pa;)“^). (2.8) 

See also (3.2) for the first few coefficients D„. The representation (2.8) generalizes [16, Eq. (2.5)] 
where the case p = 1 was considered. 

To conclude, we list the specific examples corresponding to 1 < fc < 5: 


i^'i(p, l;p,0;a;) 
iTi(p,2;p,0;a;) 

i^i(p,3;p,0;a;) 

i^'i(p,4;p,0;a;) 

iTi(p,5;p,0;a;) 


pxe^^ 


[pxfe^ ] 


{pxfe^ ] 

[- 

{px)^e^ ] 

[- 

{pxfe^ ] 

[- 


1^1 

px } 

3(1 + p) 

px 

6(1+ P) 
px 

10(1+P) 
px 


(1 + P)(2 + p) ) 

{pxY I ’ 

(1 + p)(ll + 7p) (1 + p)(2 + p)(3 + p) ) 

(px)2 + (px)3 / ’ 

5(l + p)(7 + 5p) 5(l + p)(2 + p)(5 + 3p) 

(px)^ {pxY 

(1 + p)(2 + p)(3 + p)(4 + p) 
{pxY 


3. Asymptotics of i'I'i(p, fc; p, 0; x) for fixed k 

3.1 Asymptotics for x —>■ ±oo when k is a positive integer 

When fc is a fixed positive integer and p ^ 0, the asymptotic expansion of i4>i(p, fc; p, 0; a;) as 
a: —>■ ±oo follows immediately from (2.8) as 


i4'i(p,A:;p,0;x) 


fe-i 

{pxfe^'^Dnipx)-'^ 

n—0 


(fcGN), 


(3.1) 


where the coefficients are defined in (2.7). 
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To determine the first few coefficients we note the following evaluations of the Stirling 
numbers of the hrst and second kinds: 

Then some routine algebra yields that 

Do=i, D,={1+p)D2=yi+p)(yyj{3kii+p)-i-5p}, 

D 3 = -{1 + p) {k{l + p) + 2 p}{k{l + p) + 1 + 3p},... . (3.2) 

We make three remarks concerning formula (3.1). First, it is an exact result. Secondly, this 
expansion consists of a finite series in inverse powers of x valid for a fixed admissible p 0 (that 
is, positive or negative p). And thirdly, it holds for both x —>■ +00 and x —>■ — 00 . Consequently, 
it follows that i'I'i(/9, k; p, 0; x) grows exponentially for x > 0 and decays exponentially for x < 0. 


3.2 Asymptotics for p —>■ +00 when k is a positive integer 

An approximation for i5'i(p,/c; p, 0; x) for p —>■ +00 when fc is a fixed positive integer and x is 
fixed (positive or negative) also follows directly from (2.7) and (2.8). We have 


i^'i(p,fc;p,0;x) 




k—1 n 

: ^k-n 


n=0 


e=o 





o(fe—n) 
^k-l 


p^e^ ^ 

n=0 ^ ^ 



(3.3) 


as p —>■ + 00 . More terms can be easily generated if required. 

3.3 Asymptotics for x —>■ ±00 when k > 0 

When the restriction on k being a positive integer is relaxed to A: > 0, the expansion in (3.1) 
ceases to be valid. In this case, we can apply the asymptotic theory of integral functions of 
hypergeometric type; see, for example, [2], [18, §2.3]. For the Wright function in (1.1) with the 
numerator parameters a, a and denominator parameters j3, b, we have the following (formal) 
exponential and algebraic asymptotic expansions Ei^i(z) and Hi_i(z) given by 


Ei^i(z) := Z = Kihzfl'^ (3.4) 

j=o 


and 




r((s + a)/a) 

OL ^ s! V{b — (5{s + a)/a) 


(3.5) 


where the quantities k, h and lA are defined in (1.2); see [18, pp. 56, 57]. The coefficients Aj in 
the exponential expansion (3.4) appear in the inverse factorial expansion of the ratio of gamma 
functions r(l — b + /3s)/r(l — a + as) for large |s|, with 




an algorithm for their determination is described in [18, pp. 46-48]. Then, if 0 < k < 2, we have 


i4'i(a, a;/3;&; z) 


Ei^iiz) + i7i,i(zeT-) (I arg z| < ^ttk) 
HiAze^^^) (I arg(-z)| < ^ 7 t {2 - k)) 


(3.6) 
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as \z\ —>■ c», where the upper or lower sign in is chosen according as z lies in the 

upper or lower half-plane, respectively. A more precise result that takes into account the Stokes 
phenomenon on the rays arg z = Attk when 0 < k < 1 is given in [14]. 

We now specialize the result in (3.6) to derive the asymptotics of fc; p, 0; x) for large |a;| 

in the case where p > 0 and k > 0. We have 


id'i(p,fc;p,0;a;) = ^ 

n—0 


T{k + pn) x” 
r(pn) n! ’ 


so that K = h = 1, d = k, Aq = . Then we obtain^ from (3.4)-(3.6) that 


i4'i(p,fc;p,0;a;) ~ < 


Ajx-^ 

T{is + k)/p) 
P^o r(-s-fc) 


(^+k)/P + Y Ajx-^ 
j=o 


{x —>■ -koo) 


{x —>■ —oo) 


(3.7) 

for fixed fc > 0, where the coefficients Aj = p^ ^Dj with the Dj defined in (2.7). 

The expansion of i'I'i(p, fc; p, 0; x) is seen to be dominated by the algebraic expansion as 
X —>■ —oo. We observe that this last expansion vanishes when A: is a positive integer and note the 
fact that, in this case, the coefficients Aj = 0 for j > k. The expansion in (3.7) as a: —>■ ±oo then 
agrees with that in (3.1). 


4. Asymptotics of i4>i(p, k; p, 0; x) for k —>■ -too 

In this section, we consider the asymptotics of i'I'i(p, k; p, 0; x) for k —>■ -|-oo when x is finite (or 
such that x = o{k)) and also when x is allowed to become large, namely of order k. The cases 
p > 0 and p G (—1,0) are considered separately as their treatments are distinct. Apart from 
the situation corresponding to p G (—1,0) with x finite, the approach we employ is essentially 
the saddle-point method. This is consistent with the fact that the case of negative values of 
p corresponds to the second polar type of the mechanism of formation of the probabilities of 
large deviations, whereas it is their Cramer’s or first type which is regarded as the probabilistic 
counterpart of the saddle-point method; see [22, Introduction] for more details. 


4.1 The case p > 0 

An integral representation suitable for the determination of the large-fc asymptotics of the function 
i4'i(p, fc; p, 0; a;) when p > 0 can be obtained as follows. From the loop integral representation 
for the beta function [11, Eq. (5.12.10)], we see that 


r{k + pn) _ 1 tk+pn-l 

k\T{pn) 2711 Jq {t — 1)^+1 


(k + pn > 0), 


where k > 0 (not necessarily an integer) and the integration path is a loop that starts at the 
origin, encircles the point t = 1 in the positive sense and returns to the origin. Then, substitution 
of this result into (1.3) followed by reversal of the order of summation and integration yields the 
integral representation 


i4'i(p,fc;p,0;a;) 


r(fc + pn) x” 
r(pn) n! 

n=l ' 


ki 

27ri Jo (t - 1)'”+^ 


kl /•(!+) 

2777 Jq (t - !)'=+! 


1) dt 


(4.1) 


since ^/(t — 1)*+^ has zero residue at t = 1. 

^The subdominant algebraic expansion undergoes a Stokes phenomenon on the positive real x-axis and is 
omitted. 
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4.1.1 Asymptotics for k —>■ +oo and x —>■ +oo 

To enable an estimation of k;p,0;x) for k —>■ +oo with x = 0{k), we put x = ku, where 

M > 0 is hxed. In view of (4.1), we have 


k\ f 

i'^i{p,k;p,0;x) = ^ J 


( 1 +) 


dt. 


h tit - 1) 

The function ip{t) has saddle points ts where if'it) 


ipit) = log + utP. 
0 ; that is when 


tlits - 1) 


1 

pu 


(4.2) 


(4.3) 


When 0 < p < 1, there is a single saddle^ at tso on the positive real t-axis satisfying tso > 1. 
The steepest descent path through this saddle crosses the real axis in a perpendicular direction 
and starts and finishes at the origin; see Fig. 1(a). In this case, the integration path can be 
made to coincide with the steepest descent path through tso- When p = 1, an additional saddle 
moves off an adjacent Riemann sheet and appears on the negative t-axis. When p > 1, additional 
saddles arise, which are situated symmetrically about the real t-axis on an approximately circular 
path surrounding the interval [0,1]. The directions of the valleys of ifit) (where ifit) —>■ —oo as 
|t| —>■ oo) are easily seen to be the rays arg t = ±(2s+l)7r/p, s = 0,1,...; the directions of the hills 
(or ridges) are symmetrically positioned between consecutive valleys on the rays arg t = ±27rs/p. 
Paths of steepest descent can only pass to infinity down one of the valleys or terminate at t = 0, 
whereas steepest ascent paths pass to infinity on a ridge or terminate at the singular point t = \. 

A study of the steepest descent paths through these saddles reveals that the integration 
path in (4.2) can be deformed to pass over those three saddles with the greatest real parts as 
shown in Fig. 1. The dominant saddle tso > 1 is on the positive real t-axis, with the two 
neighboring complex saddles being subdominant, since it is reasonably easy to show that 
JR(i/’(ts.±i)) < V'(^so) (we omit these details). 

Straightforward differentiation shows that at a saddle ts^ 


'^"{ts) 


2ts-l 

tiits-iy 


+ p{p - 1) utP ^ 


(1+p) (. _ 

tuts - If V i + p) 


upon use of (4.3) to simplify the second term. Then a routine application of the saddle-point 
method applied to the dominant saddle tgo (we ignore the subdominant contributions present 
when p > 1) produces 


i4'i(p, fc;p, 0;fcu) 


k\ / 2 tt exp [kipitso)] 

2Tr\l kif'itso) tsoitso - 1 ) 


{i + o(fc-)} 


= m 


27r(l -I- p) 


tsO — 


P 

l + p 


\-i/2 

j exp[k'ipitso)]{l + Oik~^)} 


as fc —>■ oo, where from (4.2) 


(4.4) 


fit so) 


log (1 - Ts) + -r 

PA - Ts) 


Ts — 1/tsO- 


(4.5) 


The value of the dominant saddle tso = tgoiu) > 1 is determined from (4.3). It is clear from (4.5) 
that fitso) > 0, with the result that the factor exp [kfitso)] is exponentially large as A: —>■ -too. 

In Table I we present some numerical values for the absolute relative error in the computation 
of iTi = i'i>iip,k; p,0]ku)/Tik) using the leading asymptotic behavior in (4.4) for different p 
and u and varying values of k. 


^It is worth mentioning that an application of [3, Proposition 5] implies (with some effort) that given p > 0, the 
solution tgo{u) to (4.3) admits the following closed-form representation in terms of the ‘reduced’ Wright function 


tsoiu) = < u - log 


e-^y 
2/(1 + 2 /) 


0-_ ^ + ^ 

p-pl yP/{p+l) 


dy 


l/(p+l) 
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Figure 1: Typical paths of steepest descent and ascent through the saddles when (a) p = 2/3, (b) p = 2, 
(c) p = 3 and (d) p = 3.6. The saddles are denoted by heavy dots and the valleys by V; the arrows 
indicate the direction of integration in (4.2) taken along steepest descent paths. There is a branch cut 
along (—oo, 1]. 


Table 1: Values of the absolute relative error in the computation of i4'i for different values of k obtained from 
(4.4). 


k 

M = 1, 

i4'i 

P=1 

Error 

u= i, 
1^1 

p = 2 

Error 

20 

1.373292 X 10^® 

1.237 X 10-2 

2.215937 X 1079 

1.084 X 10-2 

30 

1.957497 X 1027 

8.192 X 10-® 

1.200853 X 1029 

7.180 X 10-® 

50 

8.405994 X 10^® 

4.889 X 10-® 

3.021946 X 10^® 

4.286 X 10-® 

80 

6.720337 X lO^^ 

3.047 X 10-® 

3.280812 X 1077 

2.671 X 10-® 

100 

1.009961 X 10®7 

3.435 X 10-® 

7.133294 X 1096 

2.135 X 10-® 


4.1.2 Asymptotics for k —>■ +c» and k/x —>■ +cx3 

When k —+c» and x = o(fc) (which includes the case of finite x = 0(1)), we find from (4.2) or 
(4.3) the equation for the saddle points in the form 

t?(4-i) = A. 

px 

The dominant saddle on the positive real t-axis consequently has the approximate value 

tso ~ (fc-j>+oo). 

More precisely, a routine perturbation expansion procedure shows that 

’ A6) 

r=0 ' 
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Table 2: Values of the absolute relative error in the computation of I'l'i for different values of k obtained from 
(4.7) and (4.8). 


k 

X = 1, 
iTi 

P=\ 

Error 

X = 1, 
1^1 

P=1 

Error 

20 

6.966593 x 10°i 

1.117 X 10-1 

3.665626 x 10°3 

3.207 X 10-3 

50 

5.142250 X 10°2 

7.559 X 10-2 

3.555739 x lOO® 

6.793 X 10-3 

100 

3.547863 x 10°3 

5.714 X 10-2 

2.108775 X lOii 

7.198 X 10-3 

200 

3.909080 X lOO-* 

4.364 X 10-2 

3.015835 X 1021 

6.738 X 10-3 

500 

2.368065 X 10°® 

3.093 X 10-2 

5.675040 X 1036 

5.625 X 10-3 


X = 10, 

II 

X = 10, 

P=1 

k 

iTi 

Error 

i'l'i 

Error 

50 

3.511690 X lO^-* 

5.279 X 10-2 

2.620759 X IO 21 ' 

1.794 X 10-1 

100 

1.403708 X 10^3 

3.454 X 10-2 

3.486374 x 1039 

1.307 X 10-1 

200 

5.198767 X 1022 

2.311 X 10-2 

5.286625 x lO^i' 

9.396 X 10-2 

500 

4.579438 x lO^o 

1.399 X 10-2 

3.054764 X 10®® 

6.012 X 10-2 

1000 

3.340732 X lO^® 

9.751 X 10-3 

2.744393 x lOi^® 

4.273 X 10-2 


where 


-Bo — 1) Bi — 


1 + p’ 


Bo = 


2 ( 1 + 


B. = 


P(P- 1) 
3(l + p)3’ 


B, = 


p{2 — bp + 2p^) 


8(1 + P)4 


B. = 


15(1+ p)3 


(-3 + 13p- 13p2 + 3p3), 


Ba = 


144(1+ p)6 

Substitution of (4.6) into (4.5) expressed in the form 


(24 - 154p + 26V - 154p^ + 24^^),.... 


V'(iso) = 


— r 
sO 


r—1 


1 1 

- + - 
r p 


{tsO > 1 ) 


yields, with the help of Mathematica, that 


V’(4o) =V^x V 

p 2 6(1+p) 


2P-1 3 I (3p-l)(p-l) -4 

^ 10/1 ( „^9 ^ 


+ 


(4p-l)(3p-2)(2p-^ 5 ^l)(2p-l)(p-l)(p-2)^_6 , 

X. ' 


12(1+ p)2 

x"Vo(x-^) 


120(1+ p)3 60(1+ p)4 

for fc/x —>■ + 00 . From (4.4), we then obtain the leading behavior 

i4'i(p,fc;p,0;x) = -.{pxkPy/^^^^~^P'>'> exp[fc'i/'(tso)] {1 + 0{k~P^^^+P'>)}, 

\/2’r(l + P) 

as fc —>■ +00 when k/x ^ +oo, where 


kipitso) = {pk'^x 


- foi.P-T.V/(i+p) 


1 P 2^ 


, (2p-l) _2 , (3p-l)(p-l) _3 


(4.7) 


6(1+ P)^ 


12(1+ p)2 


X-VO(x-") . (4.8) 


In the cases p = 1 (where we note that the term involving y ^ jjj (4.8) vanishes) and p = 2, 
we consequently find from (4.6)-(4.8) the approximations 


1 ^ 1 ( 1 ,fc;l,0;x) = 


m 

2iV 


(kx)^/^ 
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X exp [2(A:x)^/^{l + + ^(x/Zc) + 0((a;/fc)^)}] {1 + 0(fc (4.9) 

and 

i^i(2,fc;2,0;x) = ffl(2fc2x)i/6 

VOTT 

X exp [|fc^/^(2a;)^/^{l + \{2x/k)^^^ + + 0{x/k)}] {1 + 0(fc“^/^)} (4.10) 

as fc —>■ +00 and k/x —>■ + 00 . 

From (2.1), the case p = 1 can also be expressed in terms of the confluent hypergeometric 
function iFi(l + fc; 2; x). The asymptotic behavior of this function follows with some effort from 
[20, Eq. (3.8.3)] as 

iFi(l + k] 2; x) = —— exp [2(fcx)^/^ + ^x] {1 + ^(x^/k)^^^ + 0 {x^/k)} 

as Zc —>■ +00 with k/x —>■ +00 (see also [16, Theorem 2.3] where a more accurate asymptotic 
expansion than this formula is stated for fixed x > 0). Combination of the above result with (2.1) 
then yields (4.9), albeit with a different order term. A related result on the double asymptotics of 
the function iFi can also be found in [16, Theorem 2.3.i]. Note that [16, Corollary 2.4] addresses 
the case where p = 1, Zc —>■ +00 and x = 0{k), which is not considered in the present paper. 

In Table 2 we show values of the absolute relative error in the computation of the normalized 
function I'l'i = iil'i(p, Zc; p, 0; x)/r(Zc) using the leading asymptotic behavior in (4.7) and (4.8) 
for different p and x and varying values of k. 


4.2 The case p G (—1,0) 

Throughout this subsection we consider p G (—1,0) and write p = —cr with 0 < ct < 1. From 
(3.1) and (3.2), the behavior of i4'i(— ct, Zc; ct, 0; —x) is given by 

i4'i(-ct, Zc; ct, 0; -x) ~ (CTx)''e“''|l + ^ + 

+ '"^2%xr^’' ^^{3fc(l-<^)+5a-l} + ...| (4.11) 

for X —>■ +00 and fixed positive integer Zc. From the structure of the expansion in (3.1) and the 
closed-form expressions (3.2) for the initial coefficients I?„, it follows that (4.11) with the first 
three terms in curly braces is exact for Zc = 1, 2 and 3. 

4.2.1 Asymptotics for k —>■ -|-oo and finite x 

We have for Zc > 0 (not necessarily an integer) and p = —a that 


00 

i4'i(p,Zc;p,0;x) = ^ 

n—1 


r(Zc — an) 
r(pn) 


n! 


xr(Zc) fr(Zc —ct) a^r(p) r(Zc —2 ct) x^r(p) r(Zc —3 ct) 1 
= r(fc) +2!f^ Tik) +3!f(^ r(fc) +"i- 

As Zc —>■ + 00 , the ratio r(Zc — ma)/r(k — (m + 1)ct) ~ k‘^ for finite positive integer m, with the 
consequence that 

r(Zc — ct) ^ r(Zc — 2 ct) ^ r(Zc — 3 ct) .... 

Making use of the result [11, p. 141]: 


r(Zc - ct) 

m 



cr{a + l) 
2k 



(Zc —>■ -l-oo). 


we therefore obtain^ for fixed x 

xZc^r(Zc) 


i4'i(p,Zc;p,0;x) = 


r(p) 


1 - 


P(1 - P) 
2 k 




0(fc-^)) 


®We use the shorthand notation 0(A, B, C) to represent 0{A) + 0{B) + 0{C). 
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xkPT{k) f xkPT{p) 


p(l - p) 
2 k 


^ 6r(3p) 


(l + 0(fc-i)) + 0 


x2fc2pr(p) , „ _2 \ ] 

6r(3p) ’r(2p)’r(4p)J/ 


\ 

r^; 

(4.12) 


as fc —>■ +00 for fixed x ^0 and p € (—1, 0). 

It is found numerically that (4.12) remains a good approximation when k —>■ +oo and —>■ 

+ 00 . In the case of negative x, this numerical observation can also be derived from a local large 
deviation limit theorem for partial sums of independent integer-valued random variables in the 
case of power tails of the common probability function of these variables, which is given in [1, 
Theorem 3.7.1]. The leading term in the above approximation, namely 


i'^iip,k]p,0;x) 


xk^T^k) 

~wr 


[k —>■ -foo, p G (—1,0)), 


(4.13) 


was given previously without proof by Panjer and Willmot [13, Ex. 10.2.2, p. 342]. We remark 
that when p = the refining term involving k^ inside the curly braces of the expansion in 
(4.12) vanishes; similarly, the order terms involving fc^“^/r(2p) and fc^'’/r(4p) also vanish to 
yield that 


i'I'i(-i, fc;-i, 0; a;) —--^--^=|l-I- ^[k ^+oo). (4.14) 

The leading term in (4.14) can be seen to agree with that obtained from (2.3) and (2.4) upon use 
of the large-fc asymptotic forms of the Bessel functions given in [11, Section 10.41(i)]. 


4.2.2 Asymptotics for integer k —>■ -|-oo and x —>■ —oo 

In this section, we replace negative x by —x, where x > 0. We derive an integral representation 
for i'I'i(p, fc; p, 0; —x) analogous to (4.1) when p = —a, 0 < a < 1 and fc is a positive integer. 
Now, for fc = 0,1, 2,..., we have 

r(fc -an) _ 1 ^k-an-l 

k\r{—an) 27ri J ft — 1 )*+^ ’ 


where integration is taken around a positively orientated closed path surrounding the point t = 1, 
but not enclosing the origin. This result is easily established by evaluating the residue of the 
integrand at t = 1. Then, from (1.3), it follows that 


kl 

l•(l+) 

tk-1 

2711 J 

1 

(t-i)fc+y 

fc! 

pl+) 


2711 J 

1 

(t- 


E 

n=0 

dt, 


{-xt-<^y 


dt 


n\ 


(4.15) 


where the integration path is as above. 

Suppose that x = 0(k), and set k := ux, where m > 0 is fixed. Then (4.15) can be rewritten 


as 

^1 r(i+) gkip{t) ^ 

i'i>i{-a,k;-a,0;-x) = — J V'W = log 

The phase function 'tpft) has a saddle point t^ on the positive real t-axis given by the unique root 
of the equation 

- 1) = —, (4.16) 

au 

which lies in (l,oo).^ The path of steepest descent through tg crosses the real axis in a perpen¬ 
dicular direction and terminates at the origin in a similar manner to that shown in Fig. 1(a). 


^With t = 1/y, a = au > 0, (4.16) can be expressed in the form = a(l — y) when p £ (—1,0). Then, 

by [3, Proposition 5], the unique root yg on (0,1) admits the following closed-form representation in terms of the 
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Following the same procedure as described in Section 4.1.1 (we omit the details), we therefore 
obtain that 

i4'i(-CT,fc;-CT,0;-a;) ~ exp[kip{ts)] {1 + 0{k~^)} (4.17) 

as integer k —>■ +oo and x = 0{k), where 

V'(^so) = - l0g(l - Ts)--r, Ts = l/ts. (4.18) 

'^s) 

The value of the saddle = ts{u) > 1 is determined from (4.16). 

It is seen that the leading term in (4.17) is of the same form as that in (4.4) when p > 0 (with 
p replaced by — cr), albeit with a different determination of the saddle point tg- In addition, it is 
not difficult to establish from (4.18) that 'ip{tso) < 0, with the result that the factor exp [k'ilj{tso)] 
is exponentially small as /c —>■ +oo. 


5. Asymptotics of (j){p,Q\x) for x —>■ ±oo 

In [27, 28], Wright studied the asymptotic expansion of what we have called in this paper the 
‘reduced’ Wright function (also known as a generalized Bessel function) o'I'i (—; P, z) = (j){p, S; z) 
for \z\ —>■ oo when the parameter p > 0 and p G (—1, 0). We recall from (1.4) that the particular 
case of this function corresponding to 5 = 0 is given by 

OO 

o4'i(—;p,0;a;) = (()(p,0;x) = ^ ——r, (5.1) 

nil (pn) 

n—l ^ ' 


where the variable x is real. 

In this section, we present a summary of the known asymptotics of this function as a: —>■ ±oo 
when p G (—1, 0) U (0, oo); see also [25, Appendix A.2]. It is found that when p > 0, the behavior 
of 4 >{p,0]x) is exponentially large as a; —>■ +oo, but can be exponentially large, oscillatory or 
exponentially decaying as a; —>■ —oo according to the value of p. When p G (—1,0) and x —>■ +oo, 
the behavior depends sensitively on the value of the parameter p, whereas when x —>■ —oo the 
behavior is exponentially decaying for all values of p G (—1,0). 

5.1 The expansion of (j){p,0; x) as x ^ ±oo when p > 0 

The expansion of (j){p, 0; x) when p > 0 is obtained from [27] in the form 

{ E{x) {x —>■ +oo) 

(5.2) 

E{xe'^'') + E{xe~^'^) (a; —>■ —oo). 


where 

E{x) := exp [K{hxY^^] Cj{p){K{hx)^^'^)~^ (5.3) 

V j=o 

and, from (1.2) with a = 0, /3 = p, we have k = 1 + p and h = p“^; see also [18, §2.3]. The 
coefficients Cj{p) (with co(p) = 1) are discussed in [25, Appendix A.3]. If we set 


iip) = 


(2p + l)(p + 2) 
23i pj jl 


pj{p) (j > 1), 


‘reduced’ Wright function (p for p E (—2,0): 

yg = log 


-ay „ _ 1 

</'(-(! + P). 0; 7— 


p(i - y) 


(aj/)c+l ■ 


, l/(p+l) 


(5.4) 


When p = — i, it is straightforward to verify that ys = (2a/(1 + \/Y+~4aP))^, which is consistent with the above 
formula. We do not yet have a proof when p S (—1, — although we have been able to confirm this representation 
numerically. 
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where pj{p) is a polynomial in p of degree 2(j — 1), the first few polynomials are given by 

Pi{p) = -\, p2(/i) = ^(2 - 19p +2p2), 

pjp) = —(556 + 1628p - 9093p2 + 1628^^ + 

135 

P4(p) = _ ^ (4568 - 226668p - 465702^^ + 2013479^^ - 465702^"^ - 226668p'^ + 4568p®). 

The polynomials pj{p) are related to the so-called (V.M.) Zolotarev polynomials; see [24, pp. 456- 
457] and [25, Appendix A.3]. 

When p > 1, the expansion of 4>{p., 0; x) also includes additional exponential expansions of the 
same type as E{x), but with the argument of x rotated by multiples of 27r. However, as these 
additional expansions are subdominant compared to E{x) on the positive x-axis, we can still say 
that the dominant expansion is as stated in (5.2). 

From (5.2), it can be seen that (j){p,0;x) is exponentially large as x ^ -|-oo for fixed p > 0. 
As X —>■ —oo, the growth of (/)(p, 0; x) is controlled by the exponential factors exp[K(/ix)^/'*e^’^*/'‘] 
appearing in E{xe^'^^), which depend on the sign of cos(7r//t) = cos(7r/(l -|- p)). Thus, when 
0 < p < 1 we have cos (tt/k) < 0 and (j){p, 0; x) decays exponentially as x —>■ —oo. When p = 1, 
cos(7r/«;) = 0 and (/)(p, 0;x) is oscillatory with an amplitude growing like x^!'^. Finally, when 
p > 1 we have cos (tt/k) > 0 and ())(p, 0; x) is consequently exponentially large as x —>■ —oo. 

5.2 The expansion o/())(p, 0;x) as x ^ ±oo when p G (—1,0) 

To discuss the case p G (—1,0) we set p = —cr, where 0 < ct < 1. Then use of the reflection 
formula for the gamma function shows that 

1 Tl 1 

0(—fj, 0; x) =-^ r(l + an) sinTran =-(5.5) 

TT n\ zTri 

n—l 

where 

F( 2 ;) = i4'o((t,1;—;z) ^ —r(l -kern) (0 < cr < 1). (5.6) 

nl 

n—O 

The asymptotic expansion of 4>{—(t, 0; x) as x —>■ ±oo can then be obtained from that of E{z) for 
l^l —>■ oo, which we describe below following [2], [18, §2.3]; see also [14, §3]. 

From (1.2), we have the parameters associated with the function F{z) in (5.6) given by 
K = 1 — cr and h = . We introduce the formal algebraic and exponential expansions Hip{z) 

and Ei,o{z), respectively, by 


^1.0(2) 



^-(fc + l)/<T 


(5.7) 


and 

OO 

Epo{z) := Z = ^(/iz)^/'^, (5.8) 

1=0 

where Ag = (2’kuYI'^ j n and the coefficients Cj{—a) are defined in (5.4). Then, from [2], [14], [18, 
§2.3], we have the asymptotic expansion 

( Eifi{z) + E[ifi(ze^'^^) in |argz|<7rK 

F{z) ~ { (5.9) 

[ Hip{ze^^^) in | arg(-z)| < 7r(l - k) 


as \z\ —>■ 00 , where the upper or lower sign is chosen according as z is situated in the upper or 
lower half-plane, respectively. 

In the sector | arg z\ < ^ttk, Eip{z) is exponentially large and makes the dominant contribut¬ 
ion® to F{z); see Fig. 2. The expansion Eifi{z), which is oscillatory on the rays arg z = ±i7rK 


®The subdominant algebraic expansion undergoes a Stokes phenomenon on arg z = 0. 
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Figure 2: The exponentially large (EL), exponentially small (ES) and algebraic (A) sectors associated with F(z) 
in the complex 2 :-plane when 0 < k < 1 and 9 = arg z. On the rays arg z = ztlrTK the algebraic and exponential 
(oscillatory) expansions are of comparable significance. 


(the anti-Stokes lines), is still present beyond the sector | arg z\ < ^ttk where it becomes sub¬ 
dominant in the sectors ^ttk < \ arg z\ < ttk. The rays arg z = Tttk, where Eifi{z) is maximally 
subdominant with respect to are called Stokes lines. As these rays are crossed (in the 

sense of increasing | arg z\) the exponential expansion switches off according to the now familiar 
error-function smoothing law (see [11, p. 67] for an overview of this topic) to leave the algebraic 
expansion in the sector j arg(—z)| < 7r(l — k) = ttct. The expansion on the rays arg z = ±ttk has 
been discussed in [14, §5]; see also below. 


5.2.1 The expansion o/^(—cr, 0; a;) as x ^ -boo 

From (5.5) and (5.7), the algebraic expansion associated with (j){—a,0;x) is given by 


Hix) = —{Hifi{xe 
Zm 


— 7rt(T ^TTZ 




)} 


1 “ fk+l\ Tr{k + 1) 

— / - Ti - r - sm- 

TTCT k\ \ a J a 


1 “ ^-(fe+i)/- 


(5.10) 


From (5.5) 
(j>{-a, 0 ;x) 


and (5.8), some routine algebra shows that the exponential expansion associated with 
is 


^ix) = 7^{Eifi{xe~ 

ZTTl 


^) - E^xe^^n} 


_ j^l/2^X cos(7Ta/ k) 

~ 7T ^ 


TTCT 


Cj{—a)X ^ cos < X sin-^— I - 



(5.11) 


where 

A := n{hxf/'^ (5.12) 

and we recall that k=1 — a, h = a'^. 

Let us denote the points z = appearing as the arguments of the function F{z) in (5.5) 

by P±. Then, when 0 < cr < P± lie in the exponentially large sector of Fig. 2, and consequently 
as a; —>■ -boo we have 

(j){—a,0;x)^E{x)+H{x) (0 < cr < i). (5.13) 

Since cos(7rCT/K) > 0 for this range of a, the asymptotic series E{x) in (5.11) is exponentially 
large as a; —>■ -boo and is the dominant component of the expansion. When cr = 1, P± lie on the 
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anti-Stokes lines arg z = ±^ 71 ^. On these rays, cosijra/n) = 0 and the algebraic expansion in 
(5.10) satisfies H{x) = 0, so that 




A' °° 

-7 cos 

i=o 



™3/4 °° 

^ 1=0 



2x3/2 

3V3 




(5.14) 


as X —>■ - 1 - 00 . This case is related to the Airy Ai function (see (6.3)) and, from [11, p. 198], the 
coefficients Cj{—^) can be expressed in closed form as 


A-k) 


m + i) 

54^j!r(j + i)' 


When i < (T < P± lie in the exponentially small and algebraic sectors of Fig. 2 and 

^(—cr, 0; x) ^ i?(x) -I- iJ(x) (i < (J < i). (5.15) 

Since cos {na/K) < 0 for this range of a, E{x) is exponentially small as x —>■ -|-oo and is the 
subdominant component of the expansion. 

When cr = 5 , P± lie on the Stokes lines arg z = In this case, we have the exact result: 





(x G R). 


(5.16) 


We note that cos(7r(j/«:) = —1, H{x) = 0 and, from (5.4), the coefficients Cj{—a) = 0 (j > 1) 
when cr = i. On the Stokes lines arg z = Attk, the exponential expansion of F{z) is equal to 
^Ei^q{z) together with the addition of another exponentially small expansion (see below) which 
cancels in the linear combination (5.5). This yields the result 


-E(x) =-^A^/^e ^ = ())(-|,0;-x), (5.17) 

since X = x^/d when a = ^. The fact that in this case, the ‘reduced’ Wright function coincides 
with the leading term of the corresponding exponentially small expansion is consistent with the 
well-known result that the saddle-point approximation for the probability density function of an 
inverse Gaussian distribution is exact, coinciding with this density (compare to [7, Theorem 5.6]). 
Finally, when 4 < cr < 1, P± lie in the algebraic sector and accordingly 

0(—cr, 0; x) ^ i7(x) (i < (j < 1) (5.18) 

as X —>■ -1-00. The results in (5.13)-(5.18) then describe the five different asymptotic forms of 
</)(—cr, 0; x) as X —>■ -|-oo, when cr = —p G (0,1), and agree with those given by Wright in [28]. 


5.2.2 The expansion 0 /</>(—ct, 0; x) as x ^ —00 

To deal with the expansion for x —>■ — 00 , we replace x by —x to find from (5.5) that 

(/)(-a, 0, -x) = - F(xe-"'‘)}. (5.19) 

The arguments are situated on the Stokes lines arg z = ±ttk, respectively, for all values 

of a satisfying 0 < cr < 1. It is readily seen that 

TT rri'i tj —rrifc+Ti’i', — n 

^i,o(xe ) — Jiip(xe ) = o, 

so that the algebraic expansions present in the combination (5.19) cancel. It has been shown 
in [14, Eq. (4.24)] that when the algebraic expansion iJi^o(-ze^’"’*) is optimally truncated at, or 
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near, its least term, then the exponential expansion associated with F{z) on the Stokes lines 
arg z = ±TTK assumes the form 

1 P? °° 

* ^ i=o 

Here X is defined in (5.12) and the coefficients Bj (which we do not specify here) depend on the 
coefficients Cj. It can be shown (compare [14, Section 5]) that the series involving the coefficients 
Bj cancel in the combination in (5.5), and hence, we deduce that 

0(-a,O;-x) - -^{Ei^oixe^n - 

47rz 

A' 

= ^xE^e-^J2cP-a)X-^ (0<a<l) (5.20) 

as X ^ + 00 . This expansion is seen to be exponentially small for large a; > 0 when 0 < ct < 1 and 
was obtained recently in [25, Eq. (A.14)]. The leading term in (5.20) has been derived previously 
by Mikusihski [10] who employed different methods. 

It is also relevant that a somewhat related discussion in a probabilistic context can be found 
in [29, Sections 2.10-2.11]. For instance, the leading behavior on the right-hand side of (5.14) 
can be derived with some effort from a combination of (6.3) with [29, p. 171]. 


6. Some special representations and properties of (j){p, 0; x) 


In this final section, we list some integral representations and structural properties of the ‘reduced’ 
Wright function (or generalized Bessel function) (j){p, 0; x) that we consider to be of interest. Most 
of the results we present relate to the case when p = —a, 0 < a < 1. 

6.1 Special representations 

We first note the special values of a for which 0; ±x) can be expressed in terms of a well- 

known special function. These are: 

(/)( 1 , 0 ; x) = x^/^/i( 2 a;^/^), ^( 1 , 0 ; —x) = —Ji( 2 x^/^) (x > 0 ), ( 6 . 1 ) 


()i(-i,0,x) =(x G R), 
ZyJ TT 


( 6 . 2 ) 


<(>(- 5 , 0 ,-x) 


3-^/3xAi(3-i/3x) 

_(_^r (2 {-xf/^ \ 

3V3 1 H 3V3 ) 


J—i 

3 


. 37r 


ATi 

3 


2x3/2 

3^3 


(x > 0 ), 


(^)} 


(x < 0 ) 

(6.3) 


0(-|,O,x) = (ai>0), (6.4) 

(/)(-§, 0,-x) = (a;>0), (6.5) 

where Jy, and are the usual Bessel functions, and Ai and IEk,^ are the Airy and Whittaker 
functions, respectively. The first part of expression (6.1) was noted in [25, erratum]. The expres¬ 
sions (6.3), corresponding to x > 0, (6.5) and (2.4) and the integral representations (2.5), (2.6) 
are seen to imply with some effort special cases of the structural property (1.5) for negative p. 
Similarly, the expression (6.2), when combined with the integral representation for the K Bessel 
function in [11, Eq. (10.32.10)], is also an example of (1.5). The representations (6.3) (for positive 
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x) and (6.5) have been discussed more fully in [25, Eqs. (4.9), (4.3)]; the derivation of (6.3)-(6.5) 
is given in Appendix C. 

6.2 Integral representations 

We start from the following integral representation (see [21, Eq. (2)] or [27, Section 1.4]): 

(()(—cr, 0; cc) = [ e'^~'~^'^ dr (0 < cr < 1), (6.6) 

2 Tn Jc 

where (7 is a path in the complex r-plane that starts at infinity in the third quadrant, passes round 
the origin in the positive sense and returns to infinity in the second quadrant. The representation 
(6.6) holds for complex x satisfying —tt < arg x < n. 

If the path C is collapsed onto each side of the branch cut along the negative real axis together 
with a small circular path round the origin of vanishingly small radius, we obtain the following 
representation given in [21]: 


0(—cr, 0; ±x) = =F- / exp [—± cosTTcr] sin(t'^ sinTTcr) dt (6.7) 

Jo 

for X > 0 and 0 < ct < 1 . 

Another integral representation that can be obtained from [10, Eq. (2)] is 

4 >{-a,0 ;-x) = — - - Wcr((/?) exp IVt((/?)] d(/? (x > 0; 0 < tr < 1), (6.8) 

7^(1 - 0-) Jo 


where 


:= 


sin a(f 
simp 


<r/(l rr) — a)(p 
sin ip 


An equivalent expression in the context of positive stable probability density functions can be 
derived from [29, Eq. (2.5.10)]. This is explained in [24, Eq. (3.17)], where a slightly different 
trigonometric function than Wa- (p) and integration range are used. 

A simple deduction that can be made from this last representation concerns the positivity of 
(/)(—CT, 0; —x) on the positive real axis. Since Wcr{p) > 0 for 0 < (/? < tt , it follows trivially from 
(6.8) that 

(/)(—CT, 0; —x) >0 (x > 0; 0 < CT < 1); (6.9) 


see [10, p. 192]. It is clear that the positivity stipulated by (6.9) provides an analytical justification 
for employing this ‘reduced’ Wright function in probability theory®. See also [24, Remark 3] for 
an important relationship to stochastic modelling. 

A representation of (f>{—cr, 0; —x) involving an integral of a similar function is given in the 
following theorem, which resembles [29, Eq. (2.10.II)] where an analogous formula for the Mittag- 
LefHer function is given. 


Theorem 2. For x > 0 and 0 < a < 1, we have 


(/)(-iCT,0;-x) 



/*00 

^0 


(/)(-CT,0; 


—s) ds. 


( 6 . 10 ) 


Proof. In the integral representation (6.6), we first observe that it is always possible to choose 
the path C in the r-plane such that 5i(T'^) > 0 when 0 < ct < 1. If we substitute (6.6) for 
0 ; —s) into the right-hand side of (6.10), we obtain that 


/*00 

Jo 


0; —s) ds 


1 

27ri 





dr 


2v^ 1 

X 27rt 


< 7/2 

dr = 


2^/Tr 


X 


<;i'(-iCT,0;-x) 


combination of [24, Eqs. (3.9), (3.11), (3.15)—(3.16), Remark 2] and [29, Eqs. (2.3.2)-(2.3.3)] implies with 
some effort that cr, 0; —x) > 0 for cr E [^» 1) and x E R. 
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upon reversal of the order of integration and use of the result [11, Eq. (10.32.10)] 




which is valid when 5ft (r*^) > 0. 


□ 


It is relevant that Theorem 2 constitutes an important special case of the subsequent Theorem 
3. In addition, it implies new integral representations for two specific values of the ‘reduced’ 
Wright function; see (6.11) and (6.12). Thus, if we put ct = ^ in (6.10) and use (6.2), we obtain 
that 

(()(-i,0 ;-s)ds 

Jo 

pOO 

= ^ / ,-l/2g-.V4-.V(4.) 

471 


With the change of variable s —>■ , this yields the integral representation 

^4/3 roo 

/o 


4/3 /-oo 

<?i’(-i,0;-a;) = J exp [ia;'‘/^(T'‘+ t"^)] dr. 


( 6 . 11 ) 


The validity of (6.11) can also be derived with some effort from [29, Eq. (3.4.26)]. Similarly, if 

, we obtain tl 

3-i/3x r 


we put (T = ^ in (6.10) and use (6.3), we obtain that 




2i/r 

X 

6^ 


pOO 

. / S-l/2g-.V(4.)Ai(3-l/35)ds 




ds. 


Next, we have the Laplace transforms [21, Eq. (11)] 

pOO 

/ e~^*t^~^(j){—a,6;—at~'^)dt = z~^exp{—az‘^) (0 < ct < 1, a > 0) 

Jo 


and, for p > 0, a > 0 


pOO 

/ e~^*^t^~^(j){p,6;±at^) dt = 

Jo 


z-^ exp {±az-P) {6 ^ 0) 

exp {±az~P) — 1 ((5 = 0). 


( 6 . 12 ) 


(6.13) 


(6.14) 


The result in (6.14) is obtained by substitution of the series expansion of 6 ; ±at^) followed 
by term-by-term integration. We also have the Mellin transform [21, Eq. (21)] 




/ V(-CT,0;-te '^)dt = 

/o r(o-M) 


(p > -1). 


(6.15) 


It is relevant that the evaluation (6.15) can also be derived with some effort from the moment 
representations of stable distributions (compare to [29, Eq. (2.II.11)]). 

A combination of probabilistic and analytical arguments presented in [24, Eqs. (3.9), (3.11), 
(3.15)-(3.16), Remark 2] and [29, Eqs. (2.3.2)-(2.3.3), p. 178, Property 7b] implies that for s > 0 


/ OO 

0 ; —t) dt = exp 

-OO 


(6.16) 


when cr = [1,1); the above integral is divergent when a G (0, i) on account of the exponential 
growth of (/)(—CT,0, —t) given in (5.13). This result has been verified numerically for several 
a G [^, 1). In particular, if we put cr = ^ in (6.16) and use (6.3), then we obtain for s > 0 the 
evaluation 

/ OO 

A I 7Q-l/3J.^ Jj. _ /'^3\ 1 
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A version of (6.17) is given in [29, p. 171]. It is also relevant that the functions ct,0; —t) 

with a G (0, i) belong to the so-called class of ‘trans-stable’ functions considered in [29, p. 174]. 
It is also relevant to quote [29, p. 259, Comment 2.35] and [12, p. 255, Comment after Eq. (2)], 
which implicitly claims the validity of the Fourier-transform counterpart of (6.17). 

6.3 New multiplication properties and the reflection principle 

First, we establish three general multiplication properties that involve a ‘reduced’ Wright function 
(j) expressed in terms of an integral of a product of two f functions. In view of (6.2), the first one 
generalizes Theorem 2 above. 


Theorem 3. Let cti and U 2 he such that 0 < cti < 1 and 0 < (T 2 < 1. Then for x > 0 we have 

ds 

(()(-(Ticr 2 , 0 ;-a;) = / ^(-cti, 0;-s) ())(-(J 2 , 0;-xs '”=)— . (6.18) 

Jo s 

Proof. We substitute the integral representation (6.6) for (j){—ai, 0; —s) on the right-hand side of 
(6.18) to find that 

pOC j 1 /* / /-1^\ 

I := (j){-<Ti,0;-s)(j){-a2,0;-xs~'^'^) — = — e^( 0;)dr. 

Jo s 2tti Jc \Jo s ) 

Making use of the Laplace transform result (6.13) to evaluate the inner integral, we have 

I =-^ [ dr = (j){-aia 2 , 0 ;-x) 

2Tn Jq 

by (6.6), which concludes the proof. □ 


An expression analogous to (6.18) for the class of Mittag-Leffler functions is given in [29, 
Theorem 2.10.3]. 

Theorem 4. Let pi and p 2 be positive parameters such that either (i) pi < 1, P 2 < ^ or (ii) 
pi < 1, /92 > 1 with pip 2 < 1. Then for x > 0 we have 

r°° ds 

f{-piP 2 , 0 ;-x) = (j){pi,0;-s)<l){p2,0;-xsP^) (6.19) 

Jo s 

Before giving the proof of Theorem 4, we first observe from the asymptotics of (f>{p, 0; z) in (5.2) 
and (5.3) that the absolute value of the integrand in (6.19) as s —>■ -l-oo is controlled by the 
exponential factor 


exp 


a I cos 


1 + Pi 


gi/(i+pi) _|_ cos 


1 


(2.gP2)i/(i-ep 


( 6 . 20 ) 


where aj := (1 -I- pj)~^p~^^ {j = 1,2). It is clear that, when pi < I, P 2 < 1, or pi < 1, 
P 2 > I (with P 1 P 2 < I), this factor decreases exponentially as s —>■ -|-oo, so that the integral in 
(6.19) converges absolutely. However, when pi > 1, p 2 > I, or pi > 1, p 2 < 1 (with pip 2 < I), 
the exponential factor (6.20) increases without limit as s —>■ -l-oo, with the result that the integral 
in (6.19) is not defined in these cases. 

When Pi = p 2 = 1, the exponential factor (6.20) equals 1 and the (improper) integral on the 
right-hand side of (6.19) does not exist. This can be understood by use of the second relation in 
(6.1). The integral in (6.19) can then be expressed in terms of 

tJi{ 2 t)Ji( 2 x^^^t) dt. 

But in view of [11, Eq. (10.22.56)] this integral does not converge. 



Proof. From the representation analogous to (6.6), which can be found in [27, §1.4], one obtains 
that for p > 0 


(/)(p,0;-x) 


1 

27ri 


dr. 


( 6 . 21 ) 
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where C is the loop surrounding the negative real r-axis as before. The right-hand side of (6.19) 
can therefore be written as 

^1 = ^/ e^(^J 

The inner integral may be evaluated by use of the Laplace transform in (6.14) to yield that 

A = ^ ^ - 1} dr = ^(-pipa, 0; -x) {pip 2 < 1) 

by ( 6 . 6 ), where the —1 in the integrand makes no contribution since e'^dr = 0. The condition 
P 1 P 2 < 1) which is clearly satisfied when (i) or (ii) is fulfilled, applies provided the integral in 
(6.19) converges. This yields the conditions stated in the theorem. □ 


Theorem 5. Let a and p be such that 0 < cr < 1 and p > 0. Then for x > 0 we have 

ds 

(j){pa,0]±x) = / <p{—a,0]—s)(j){p,0]±xs ^)— . (6.22) 

Jo s 


Proof. It follows the same procedure employed in establishing Theorem 3 and uses ( 6 . 6 ), the 
Laplace transform result (6.14) and, finally, the representation (6.21). □ 

Corollary 1. For a € (0,1) and x > 0, the following reflection principle holds: 

(j){a,0]x) = \fx / (/)(—CT, 0; —s) s“^'^/i(2-yiEs) ds 

do 

> 

poo 

(/)(—CT, 0; —cc) = —-x/x / <))(ct, 0; —s) s“^/^ Ji(2-\/(es) ds. 

do 


Proof. If we put p = I in (6.22), we find that for a; > 0 

ds 

(/)(cr, 0; ±x) = / ())(—cr, 0; —s) (/)(1,0; ices)— (0 < ct < I). 

do s 


(6.23) 


(6.24) 


The result (6.23) follows immediately by taking the upper signs in (6.24) and using (6.19) with 
Pi = cr G (0,1), p 2 = I combined with the Bessel function representations in (6.1). □ 


Some other integrals involving products of ‘reduced’ Wright functions are cited in [21]. We 
remark that (6.18) can also be obtained with some effort from [29, Eq. (3.3.1)] which states^ that 
Li(cri)L 2 (cr 2 )=T(cricr 2 ), where Li(cri) and 42 (o’ 2 ) are two independent positive stable random 
variables with corresponding indices of stability ai and CT 2 , and T (< 71172 ) is the positive stable 
random variable with index of stability criCT 2 such that the corresponding p.d.f.’s of these random 
variables are given by x~^<j){—a,0-, —x~^) (compare to [24, Eq. (3.10)]). 

A combination of the representations (6.1)-(6.3) and (6.5) with <7 = 5,5 and | then leads to 
the following integral representations in terms of products of well-known special functions: 


1 Hr 

^(i, 0; ±x) = ±-y — J /"‘Ci(2-\/sx) ds, 

fx /2s^/^\ 

(j){l,0,±x) = ±—J sKly-^^jCl{2^/^)ds 

and 

</)(§, 0 ,±a;) = (^) 

^The sign = means that the distributions of random variables coincide. 


(6.25) 

(6.26) 

(6.27) 
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for a; > 0, where Ci denotes the Bessel functions /i (upper signs) and Ji (lower signs). 

We regard Corollary 1 and its illustrations given by (6.25)-(6.27) to be of interest in their own 
right. However, it is yet unclear whether these formulas might have a computational value. For 
instance, a substitution of the well-known series representations for Ii and Ji into (6.25) shows 


that 


4>{\,Q]±x) 


1 ^ r(in-I-i)(±2a;) 

2y^ r(n)n! 

'' n—1 ^ ' 


which is comparable to the original series representation of this function given by (5.1). 

Our final theorem concerns an asymptotic connection between the Wright function iiki and 
its reduced counterpart (f). 


Theorem 6. For p > 0 and x > 0, we have as k ^ -l-oo 


i'i'i{p,k;p,0;x) ~ r{k) 4>{p,0; xk^). 


(6.28) 


The same result holds when p € (—1,0) and x € R as fc —>■ -l-oo. 

Proof. Consider first the case with p > 0 and x > 0. With k = 1 + p and h = p^, the leading 
behavior of i4'i(p, k;p,0-,x) from (4.7) and (4.8) can be expressed in the following form: 

I'i’iip, k; p, 0; x) ~ r(fc) ((hxk^)^^'^'^ exp [K{hxk^)^^^] 

\ ZTTK J 

as fc —>■ - 1-00 with k/x —>■ -l-oo. But from (5.2), the leading behavior of (/)(p, 0; xk^) is 

Y (/ixF)^/2k gxp [K(/ja;fcP)i/«l 
TTK / 

as xkP —>■ -l-oo, which implies (6.28). 

When p G (—1, 0), we have from (4.13) with p = —a, 0 < a <1 and x G R that 
T / 7 n \ xfc“'^r(fc) 

i^'i(-(T, fc;-cr, 0;x) ~ ^ (fc-;>-l-oo). 

1 (-cr) 

From the series representation for 4>{—a, 0; xk~'^) in (5.1), it is readily seen that as fc —>■ -l-oo, the 
argument xk~'^ —^ 0 for fixed x, and the first term yields the leading behavior. Hence, we have 


^(-cr, 0; xk~'’) ~ —-- {k -koo) 

r(-cr) 

and the result follows. □ 

To conclude this section, we note that an asymptotic relationship of the same type as (6.28) 
in the case where p G (—1,—^], k —>■ -foo and x —>■ — oo can be derived with some effort from 
a probabilistic local limit theorem stipulated by [8, Eq. (1.8)]. (Note that the case p G (—i,0) 
should be excluded in order to have the finiteness of the so-called first pseudomoment; compare 
[22, Section 3.2].) 


Appendix A: Proof of the structural representation (1.5) 


In this appendix, we establish the structural representation given in (1.5) that relates to the 
‘reduced’ Wright function 4>- This can be separated into two cases according as p is positive or 
negative as follows: 


i^'i(p,fc;p,0;x) 


i^'i(p,fc;p,0;-x) 




e ^(j){p,Q\XT'’)dT (p > 0, k > 0) 
t^(p,0;—XT^) dr (pG(—1,0), 


fc = 0,l,2,...). 


(A.l) 
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The case when p > 0 is straightforward. Substitution of the series expansion for (j){p, 0; xt^) in 
(5.1), followed by reversal of the order of summation and integration, yields that the right-hand 
side of the first expression in (A.l) can be written in the following form: 


E 


n! r(pn) 



OO 

e-^^+P’^-idr = 

n—1 


r(fc -I- pn) a;" 
r(pn) n! 


valid for k > 0. By (1.3), this is seen to equal the left-hand side of (A.l). We note that the 
case k = 0 is trivial, since i’l'i(p,0; p,0;x) = and, from (6.13), the first integral in (A.l) also 
reduces to e^. 

When p G (—1,0) we put p = — cr in the second expression in (A.l), with 0 < ct < 1. We then 
have to show that 

poo 

id>i(—(j, fc; —cr, 0; —x) = / dr (A.2) 

Jo 

for non-negative integer k. Now, for fc = 0,1, 2,..., we have 


i'i'i{-a,k;-a,0;-x) 


•s^ r(fc — an) (—x)" 

^ r{—an) n! 

n—1 ^ ' 

/ \n 

(—)^ E —— 1)... {an — fc -I- 1), (A.3) 


with the product being empty (equal to 1) when fc = 0. But from (6.13) we have the Laplace 
transform for z > 0 

4>{—a, 0; —xt~'^) dr = exp (—xz'^). 

Compare to [24, p. 452] where a connection to a positive stable density function and its Laplace 
transform is made; see also a comment above Theorem 6. Differentiation of this result k times 
then yields that 



pOO 

Jo 


(j){—a,0]—XT '^) 


dr 


ik °° rjk 

n—1 


(-fE 


{-xf 


an{an — 1)... {an — /c -I- 1). (A.4) 


The evaluation in (A.4) is seen to equal the series representation (A.3) of the left-hand side 
of (A.2) when fc = 0,1, 2,..., thereby concluding the proof; compare to [4, Theorem la]. □ 


Appendix B: Proof of the integral representations in Theorem 1 


Denote the integral on the right-hand side of (2.5) by From the Mellin-Barnes integral repre¬ 
sentation for the modified Bessel function of the second kind (see, for example, [11, Eq. (10.32.13)]) 
we have 


Ki{z) 

3 ' ' 


4711 


^C+CXDZ 


T{s)T{s-mz)-^^ds 


(c > i, |arg z| < \tt). 


With z = ajs/t^ a = 2(x/3)^/^ (x > 0), we then find that 

7-3/2 


h = 


377 


X 


pOO 

/ e~*t'^~^^‘^Ki{a/Vi)dt 

Jo ^ 

*00 p -I pc-\-ooi N 

^-t^k-513 I — y r(s)r(s - i)(x/3)"3"C(is j dt 

-| pC-\-OOl 

— r{s)T{s - i)r(fc + s - l){x/3)-^-^ds 

6v37r 2TnJ^_^i 


3^/^7r 

x^ 
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upon reversal of the order of integration and evaluation of the integral over t as a gamma function. 
Making use of the triplication formula 


r(s)r(s-i)r(s-|) = 


27r 


-r(3s-2) 


3 VI-' 3/ 32S-5/2 

for the gamma function (see, for example, [11, Eq. (5.5.6)]), we obtain that 


(B.l) 


3a; 


^c+ooi 


h = —J r(3s-2) 

Li n V I ^ — CVl'J 


r(fc + s — |) _ 




X 


For k = 1,2,..., the integrand has poles only at s = 3 — |n (n = 1,2,...) resulting from 
r(3s — 2), since the quotient of gamma functions is regular. Displacement of the integration path 
to the left to coincide with the vertical line JR(s) = | — ^N, passing over the first iV — 1 poles of 
r(3s — 2) with residue |(—)"'/n!, we find 


_ r(fc gu) (—a;)’' 

^ ^ r(-in) n! 


Rn, 


(B.2) 


n—1 


where, on putting s = | — with —00 < t < 00 , 


Rn = 


6 3 

A-l/2 poo 


2 tt 


/ oo k — 1 

r(-V+i+lt) C/Ar(t)x“*‘(it, GN{t) = ]^(s-|+m). 

m=0 


Now, |s — I + m| < ^{N'^ + provided N > 3m > 1, so that 

\GN{t)\<{N/3)'^il + t^)'^/^ {N>3k). 

By repeated application of the well-known formula zr(z) = r(z -|- 1), we see that 


|r( — N — 


\nh+m 


|(V—i— it)(iV —| — it)... ( 5 —it)| 
With these bounds we then have 

a;^-i/2(V/3)'= 


< 


m+h) 


liivl < 


2 v^r(iv + i) 


|r(i + zt)|(i + t2)'=/2(it, 


where the integral is convergent, since |r (5 -I- it)| = (tt/ coshTrt)^/^. It therefore follows that 
Rn = 0 {N^x^~^/^)/r{N + i) and hence, that Rn —>■ 0 as V —>• 00 for fixed a; > 0. We can then 
put iV = 00 in (B.2), which is seen, by (1.3), to be the series expansion of i4'i(—i, fc; —i, 0; —a;). 
This establishes (2.5) for positive integer k. 

A similar procedure can be employed to establish (2.6). We have the Mellin-Barnes integral 
representation given in [11, Eq. (13.16.12)] 


i(z) = — [ 

V 6 2Tri 


c-\-ooi 


r(,+ i)r(, + |) __ 


z ^ds (c > 0, |argz|<i7r). 


r(s + i) 

With z = a = ^x^ 121 {x > 0), the integral on the right-hand side of (2.6) may be expressed 


as 


li = J e ^ exp(—ia/t^)IEi 1 (a/t^)dt 


3 ^ fC+ooi l)ris+l) 


TT 27ri 


C—OCl 
C+002 




r(s + i) 

r(fc-b 2 s) 
r(2s) 


r(fc-b 2 s) 




ds 


3s 


ds 
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upon use of the duplication and triplication formulas for the gamma function; compare to (B.l). 

Provided k = 1,2,..., the only poles of the above integrand result from the factor r(3s) 
situated at s = —^n {n = 1,2,...). Displacement of the integration path to the left to coincide 
with the vertical line 5ft(s) = —^N + i and evaluation of the residues then yields that 


h = 


^^r(ft-|n) {-xY 

k r(-in) 


n! 


Rn, 


(B.3) 


where, by similar arguments to those above, Rm = /T{N + i) —>■ 0 as fV —>■ oo. 

Hence, we can put = oo in (B.3), which then equals the series expansion of |, fc; —|, 0; —x). 

This establishes (2.6) for positive integer k. □ 


Appendix C: Proof of the special function representations in Section 6.1 

From (5.1) and (5.5) with a; > 0, we have 

1 / I \ ji 

0; zbx) =--;— r(l + ^n) sin(^7rn). 

n—l 


Replacing n by 3m + j, j = 0,1, 2, we obtain 

r(m + |) 
(3m + 1)! 


</)(-i,0;±a:) = T^{f; 


(Tx) 

(3m + 2)! 

m—O ^ ' 


3m 


since the sum corresponding to j = 0 vanishes. Use of the triplication formula (see (B.l)) for the 
gamma function and introduction of the new variable X = x^ jXI leads to 


(l}i-^,0;±x) = ^ 


(txy 


m—0 


m!r(m + 




{Txy 


m—0 


m\T{m + 


X 

= Ttt 


(2Vt) + A’i/®Ci(2yA) j. 


(C.l) 


where C^{z) denotes the Bessel functions Ju{z) (upper signs) and Iv{z) (lower signs). 
Taking the upper and lower signs in (C.l) separately, we therefore find that 


()i(-i,0;x) = - 


3V3 




and 




(C.3) 


The representation of (C.2) and (C.3) in terms of the Ai Airy function follows upon use of [11, 
Eqs. (9.6.1), (9.6.6)]. This completes the proof of (6.3). 

The proof of (6.4) and (6.5) follows the same procedure. From (5.5), we have 


(j}i-l,0]±x) = --Y^ 


{±x) 

nl 


■ r(l + |n) sin(|7rn). 


which becomes, on replacing n by 3m + j, j = 0,1, 2 

XyJZ 


(/)(-|,0;±a;) = T' 


27r 




m=0 


m—0 


(3m+ 2)! 
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Use of the duplication and triplication formulas for the gamma function and introduction of the 
new variable X = 127 then leads to 




= T' 


301 


= T 


= T 




T2: 


(±v)' 


^m=0 

.22/3^^ rr(|) 


3v^ lr(|) 
rr(i) 


r("i+3) 


6v^ lr(| 


f(iUi(|;|;±A’)±A’i/3 


3) ^ P / 7. 4 


r(|) 




(C.4) 


From [11, Eqs. (13.2.42), (13.14.5)], we have the following combination of confluent hyperge¬ 
ometric functions expressed in terms of the Whittaker function 1 Fk_^(z) given by 


r(i-fe) 

r(a- 6-b 1 ) 


iFi{a-,b]z) + 


r(fe-i) 

r(a) 


iFi(a 


6+1; 2 


b;z) 


= (C.5) 

for z G C. Inserting the values a = |, 6 = | in (C.5), we find upon taking the upper signs in 
(C.4) that 


0(-|,O;a:) = ei^X-^/^ W_,.{X) 


}_^2F/27w_, J 

2v^ "’"V 


2 ’ 6 
4a;3 


{x > 0), 


(C.6) 


where we have made use of the fact that IF„,_^(z) = Wk^^{z). This completes the proof of (6.4). 
Taking the lower signs in (C.4), we have 


22/3x 

</)(-|,0 ;-x) = ^e 




upon use of the well-known Kummer transformation 1+1(0; b; z) = 1+1(6 — o; 6; —z). Insertion 

of the values a = —g,6 =|in (C.4) then yields that 


(/)(-§, 0;-a;) = 


2^ 

■\/37r 


e-x/2x-i/3w^ i(A') 


= r( — ) (x>0). 

V TT V 27 y ^ ^ 


(C.7) 


This completes the proof of (6.5). □ 

Although we do not use this result here, it is worth recording the connection between the 
Whittaker functions W^i i{z). From [11, Eqs. (13.15.8), (13.18.0)], we obtain the relation 


f)2, 

IU_i 1(2) = 6IFi 1(0 )-^A:i(iz) {zGC). 
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